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Abstract 

We have discovered an unexpected and surprising fact: a 2D axially symmetric short-range 



potential contains infinite number of the levels of negative energy if one takes into account the spin- 
orbit (SO) interaction. For a shallow well (m e UoR 2 /h 2 <C 1, where m e is the effective mass, Uq and 
R are the depth and the radius of the well, correspondingly) and weak SO coupling (\a\m e R/H <C 1, 
a is the SO coupling constant) exactly one two-fold degenerate bound state exists for each value of 
the half-integer moment j = m + 1/2, and the corresponding binding energy E m extremely rapidly 
i decreases with increasing m. 
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As it is well known from any textbook on quantum mechanics a very shallow potential 
well (m e UoR 2 /h 2 <C 1) cannot capture a particle with the mass m e in 3D case and does 
this in 2D and ID situations provided the wells are symmetric: even potential in ID, axially 
symmetric well in 2D. In the latter case the only negative level corresponds to the s-state 
(m = 0). 

Consider now a 2D electron with accounting for the SO interaction in Bychkov-Rashba 
form 1]; the Hamiltonian reads: 

"2 

where r and p are the radius in cylindrical coordinates and the 2D electron momentum 
operator, respectively, <r are Pauli matrices, n is the normal to the plane of 2D system. 
It is convenient to write down the Schrodinger equation in the p-representation: 

+ a(<x[p x n])]¥(p) + J ^U(p - p')¥(p') = 5*(p). (2) 

Here U(p) = J dre~ ipr U (r) = 2n J °° drrll (r) J (pr) is the Fourier-transform of the potential 
(Jo(z) is the Bessel function). Because of the axial symmetry of the problem it is possible 
to separate the cylindrical harmonics of the spinor wave function and search for the solution 
in the form 

*M(p) = { f*<P*» ) , (3 ) 
(if is the azimuthal angle of the vector p). Using the summation theorem 

oo 

^o(|p-p'IO= MP 1 ") J o(p'r) cos (k6) 

k=—oo 

(8 is the angle between the vectors p and p') one can rewrite Eq.(J2J) for each m-th harmonic: 
1 £}4"?(P) ± iocp^ip) + / dr'r'U(r')J mHlTl)/2 (pr')Cif(r') = 0. (4) 
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Here the functions C-[ ™ have been introduced: 

POO , - 



C{ ; (r) = / dppj m (pr)ip\ \p), C\ >(r) = / dppJ m+ i(pr)^ { 2 \p). (5) 
Jo Jo 

Resolving Eq.Q with regard to ^1.2 we find: 

^\P) \ 1 r jT r^( [P 2 /2m e - £]Jm(pr)C{ m \r) - tapJ m+1 (pr)Ci m \r) 

drru [rj 



4 m \p) / A(p; \ ta P J m (pr)C[ m \r) + \p 2 /2m e - £}J m+l (pr)C { ™\r) 



where A(p;£) = (p 2 /2m e — £) 2 — 4a 2 p 2 . Zeros of A(p; £) as functions of £ give two branches 
of the dispersion relation for free electrons: £±(p) = p 2 /2m e ± ap. 

Finally, from Eq.fjUJ) using the definitions © we arrive at the equations for C["^ 

POO 

C { r\r)= dr'r'U(r')A^\r,r')d m \r') (*,J = 1,2). (7) 



Here the matrix A-"^ has been introduced: 



A^(r,r') = -ia ^—J m ipr)J m+1 (pr'), A^(r,r') = (A^(r',r)y. (9) 

The function A(p; £) can be presented in the form A(p; E) = \{jp — p ) 2 /2m e — E] ■ [(p + 
Po) 2 /2m e —E], where po = m e \a\ is the radius of the loop of extrema, E is the energy counted 
from the bottom of continuum, E = £ + m e a 2 /2. Now we search for levels of negative energy 
satisfying the condition \E\ <C m e a 2 and simultaneously we assume 2m e UoR 2 /h 2 = ^ C 1 
(Uo, R are the characteristic depth and radius of the well). Then integrals in Eqs.(|BEJ) can be 
calculated in the "pole" approximation: we put p = Po everywhere in the integrand except 
the first factor in A(p; E). As a result we have (h = 1): 
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o 



-i sign(a) Jm{p r)}J m+1 (p r')C { 2 m) (r')} 



Ct\r) = / dr'r'U{r')[i sign(a) J m+1 ( Po r)} J m (Por')Cl m) (r') 

y l /2\E\ Jo 

+J m+ i (Por)} J m+ i {pvr')C { ™ ] (r')] (10) 

Thus, we obtained the system of linear integral equations with degenerate kernels which can 
be easily solved. This system can be reduced to a pair of linear algebraic equations for the 
quantities t m = J °° drrll (r) J m (p r) and t m+ i (defined similarly): 

t m = j== [t m ~ i sign(a)t m+ i], 
V 2 I-^I 

t m+1 = ^L [t m+1 + i sign(a)t m ]. (11) 
Here Xm = —npn^/ml J °° drrU {r) J^Pqt) . From Eq. ffTT]) one immediately gets 



E m = -(Xm + X m+ i) 2 /2 = -t^l ( r dr rU(r)[J 2 m (p r) + J 2 m+1 (p r)f ' 



o 



„2 2 / /-oo \ 2 

' '■" ' ' 1 dr rU(r)[J 2 _ 1/2 (p r) + J 2 +1/2 (p r)}) , (12) 



o 



where j = m+1/2 = ±1/2, ±3/2... is the z-projection of the total moment. As it is seen from 
Eq. (fT2*|) all levels are two-fold degenerate: E m is even function of j. If now SO interaction 
is small (poR *C 1) we can get the asymptotic behavior of the binding energy by expanding 
the Bessel functions in Eq. (jl2|) . For a rectangular well U(r) = —Uq8(R — r) we have 
\E m \ oc a 4 l J 'l/(2 4|il ((|j| - l/2)!) 4 (2|j| + l) 2 ). For an exponential well U(r) = -U exp(-r/R) 
one can find: \E m \ oc a 4 W((2|j|)!) 2 /(2 4 W((|j| - 1/2)!) 4 ). 

Thus, we see that in an arbitrary axially symmetric short-range (the integral in Eq. (jl2}) 
converges) potential well there exists at least one bound state for each cylindrical harmonic 
with the energy level below the bottom of continuum (— m e a 2 /2). The energy of this state 
E counted from —m e a 2 /2 in the regime \E\ <C m e a 2 is proportional to Uq, where Uq is 
the characteristic depth of the well. Such a dependence is typical for a shallow level in 
a symmetric ID potential well. One-dimensional character of the motion results from the 
so called "loop of extrema" (see j^|). In a small vicinity of the bottom of continuum the 
dispersion law of 2D electrons has a form £(p) = —m e a 2 /2 + (p — po) 2 /2m e and corresponds 
to a ID particle at least in the sense of the density of states: one may formally consider 
the problem as the motion of a particle with anisotropic effective mass; in the p-space the 
radial component of the mass equals m e , while its azimuthal component is infinitely large 
(the dispersion law is independent of the angle in p-plane).j4| 

We realize that our conclusion looks paradoxically: for a sufficiently large value of m the 
centrifugal barrier (CB) can make the effective potential energy U(r) + Ucb positive all over 
the space. How can a bound state with negative energy be formed in such a situation? Our 
arguments are as follows: for a particle with dispersion relation (p — po) 2 /2m e there exists 
no CB; the azimuthal effective mass tends to infinity and CB vanishes. 

To check our results we have numerically analyzed the square well potential U(r) = 
—Uq6(R — r) where 9 is the Heaviside function. We seek for a solution of Schrodinger 
equation in the form 

where now if is the azimuthal angle of the vector r. Spinor components ^1,2 (r) are 
given by linear combinations of the Bessel functions J m (k±r), J m+1 (k±r) for r < R or 
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K m (K±r), K m+ i(K±r) for r > R, where 



k± = \j2m e (E + Uq) ± ma; 



k± = \/2m e \E\ ± ima. 



(14) 



Expressions (jUj) are valid when the condition E < is satisfied. Now we have to meet 
the matching conditions for the wave function and its derivative at r = R. After rather 
cumbersome algebra we arrive at the determinants, zeros of which give the required spectrum 
of localized states. The energy levels have been estimated numerically for s- and p-states 
(m=0, 1). The results totaly coincide with the ones given above for \E\ <C m e a 2 . Fig.l 
demonstrates this for s-state. The exited p-state at zero SO interaction appears when Uq 
exceeds a certain critical value Uq C \ namely, when £ > £ c = x\, where x\ is the first root 
of the Bessel function Jq(x). Taking into account SO interaction results in splitting of the 
p-level and lowering the critical value Uq for the upper of spin-split sublevels. The lower 
sublevel exists at any value of the parameter £ (see Fig. 2). 

Our last remark relates to the Coulomb interaction with a charged impurity. If one tries to 
apply the general relation (jl2j) to the Coulomb potential —e 2 /r the integral logarithmically 
diverges at the upper limit: J^z — > oo) ~ (2/irz) cos 2 (z — Tim/ A — 7r/4) does not depend 
on m after averaging over oscillations and can be replaced by 1/ttz. Eq. (jl2)) gives for the 
energy the value independent of m: 



where L is some cut-off length. Its exact value depends on the concrete situations: it may 
be the screening radius or the thickness of 2DEG. We see that the energy spectrum does 
not depend on m (as it must be for the Coulomb field) and exactly coincides with that of a 
"ID hydrogen atom": the ground state binding energy equals 2Ry (rather than Ry/2 as in 
3D case) multiplied In 2 (A), where A is the cut-off parameter (see, for example, the problem 
of the hydrogen atom in an extremely high magnetic field P). This result also supports our 
interpretation: in the region \E\ -C m e a 2 the particle becomes effectively one-dimensional. 

It is interesting from the general physics point of view to find a similar situation for 
3D case. E.G.Batyev has kindly reminded us that the roton spectrum of liquid He-4 also 
contains a part of dispersion relation that reads A + (p — p ) 2 /2M, possesses not a loop 
but a surface of extrema, and, correspondingly, should describe an effectively ID particle. 



E 



—2m e e 



4 In 2 ( Po L), 



(15) 
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FIG. 1: The behavior of the s-level versus the well depth. The curves demonstrate the transition 
between 2D and ID regimes. At a = we have an exponentially shallow level (2D result), while 
for finite a at small enough £ = 2m e UoR 2 the binding energy parabolically depends on Uq (ID 
regime). The dotted line represents the results of our "pole" approximation (Ea.p2jl) 

We have made the proper calculation, in other words we solved the Schrodinger equation 
in the momentum representation for the Hamiltonian A + (p — po) 2 /2M + U(r) with U(r) 
as an attractive spherically symmetric potential. We used the same method - expansion 
of the wave function over the spherical harmonics and we got the same result: even in 
3D a shallow potential well contains one bound state for each moment I and this state is 
(21 + l)-fold degenerate: 

E = -2n 2 p 2 M QH dr rU(r)J 2 +1/2 (p r)^ . (16) 

For the Coulomb potential the last formula once again leads to the ID result given by 
Eq.CESJ). 

In conclusion, we have shown that 2D electrons interact with impurities by a very special 
way if one takes into account SO coupling: due to the loop of extrema the system behaves 
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FIG. 2: p-states. Comparison of the exact solution for the square well with Eq.(|12|) (dotted line). 

(c) 

Inset: spin split states of the p-level: the upper curve terminates at Uq = Uq - the level merges 
with continuum. 

as ID one for negative energies close to the bottom of continuum. This results in the infinite 
number of bound states even for a short range potential. 
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